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Abstract 

We analyze an idealized model for the transmission or flow of particles, or discrete packets of 
information, in a weight bearing branching hierarchical 2 — D networks, and its variants. The 
capacities add hierarchically down the clusters. Each node can accommodate a limited number 
of packets, depending on its capacity and the packets hop from node to node, following the links 
between the nodes. The statistical properties of this system are given by the Maxwell - Boltzmann 
distribution. We obtain analytical expressions for the mean occupation numbers as functions of 
capacity, for different network topologies. The analytical results are shown to be in agreement 
with the numerical simulations. The traffic flow in these models can be represented by the site 
percolation problem. It is seen that the percolation transitions in the 2 — D model and in its variant 
lattices are continuous transitions, whereas the transition is found to be explosive (discontinuous) 
for the V— lattice, the critical case of the 2 — D lattice. We discuss the implications of our analysis. 
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I. INTRODUCTION 



Studies of dynamical behavior on complex networks have witnessed exponential growth 
in recent years. This rapid growth has been fueled by the application of network science 
to diverse areas like biology, sociology, physics and computer science and others P-[3]. The 
study of the flow of traffic through a network is a important and frequently encountered 
dynamical process. The dynamics of the flow of information has been studied over a wide 
range of communication networks of various types [il-ITT] , with the Internet being an impor- 
tant example. In the internet, computers of different computational processing capacities 
communicate by exchanging information in the form of discrete units, or "packets" which 
then find their destination by hopping from node to node via the links between them. Thus, 
the transmission of information on the network is mapped onto the traffic of packets or 
particles moving in the network via connections. The packets or particles are often assumed 
to be non-interactive, however, to understand phenomena like congestion, the interaction 
between the packets has to to be taken into account. In real networks, the finite processing 
capacity or processing power of nodes causes collective behavior such as congestion, jam- 
ming or the failure of transmission of packets. The effect of congestion in networks has been 
studied via centrality measures such as the betweenness centrality [T2HI7j, as well as via 
order parameters like the fraction of packets delivered (TUl IH [IE, El]- Cellular automata 
models have also been useful in studying the dynamics of the transmission of particles or 
packets in the network [T5ti25] . 

An interesting direction of analysis for packet transport has been to explore the stationary 
distributions of particles at various nodes, and to map them on to the usual distributions 
of statistical ensembles [26l [27] . This has been carried out both for lattices of homogeneous 
capacity at nodes and for the more realistic case of inhomogeneous capacity at nodes. The 
model [26], in which the particles move randomly in the network with the constraint that 
each node can be occupied by at most one particle, conforms to the Fermi -Dirac distribution. 
Another model [2^, in which particles performing a random walk on nodes of finite process- 
ing capacity with the constraint that no particle can move to a node which has already has 
its maximum permitted number of packets, satisfies the Maxwell-Boltzmann distribution for 
the statistics of particles. Most of the earlier studies of the transport of packets have been 
carried out for scale-free or random networks. However, very little analysis has been carried 
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out for the transport of packets in tree- or river- like branching hierarchical networks. In 
this paper, we explore the transport dynamics and statistics on a hierarchical regular 2D 
lattice model, two of its variants [2H], and its critical case [2^. The base 2D model is, in 
fact, the (7(0, 1) case of the Coppersmith model [30]. Analogous models are river network 
models [31], voter models [32j, the Domany Kinzel cellular-automata model [33] and the 
branching hierarchical model of the lung inflation process [311 ES]. In all the four lattice 
networks analyzed by us, the nodes bear finite processing capacity and the packets hop from 
node to node utilizing the lattice links, incorporating the constraint that no packet can hop 
to a node that has already has its maximum allowed number of packets. If the capacity 
distribution of a node is mapped to energy levels, the packet dynamics of transport in these 
hierarchical networks is found to follow the Maxwell - Boltzmann distribution. 

The transport of packets on a network is also analogous to the site percolation problem. 
When the percolating network supports a giant cluster, then communication is possible 
through this large connected component, otherwise communication breaks down. Recently, 
Achlioptas et al. [36] studied percolation for the Erdos Renyi model using a product rule and 
found that the giant component emerged suddenly at the percolation threshold, and that the 
percolation transition was discontinuous. This discontinuous percolation transition appears 
when the growth of the largest cluster is systematically suppressed thereby promoting the 
formation of several large components that eventually merge in an explosive way [37] . Several 
aggregation models, based on percolation, have been developed to achieve this change in 
the nature of transition. [36], |38H11] • Recent studies argue that the transition is really a 
continuous transition, with the discontinuity seen being due to the finite size of the lattice 

[laile]. 

The explosive percolation transition is seen in our hierarchical models in an unusual con- 
text. The percolation transitions in our 2D model, and its two variants, have also been 
explored and are found to be the usual continuous transition. However, the percolation 
transition in the V— lattice, the critical case of the 2D model is found to be discontinu- 
ous. This explosive transition does not follow the scaling laws of the continuous percolation 
transition. It is interesting to note that the explosive percolation transition is seen in only 
one realization, viz. the V— lattice or the critical realization of the original lattice. This 
realization is also the only case where power-law scaling is observed in the avalanche distri- 
butions. Thus the structure of the V— lattice appears to contribute to the critical behavior 
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FIG. 1. The (a) original lattice, (b) the V— lattice, (c) the strategy-I lattice and (d) the strategy-II 
lattice of M = 4 with 4 sites per layer. The numbers in brackets beside each site denote the weight 
bearing capacity of that site. The thick line in each panel denotes the trunk of the respective 
lattice (network). 

of transport processes on the structure. The analysis of this case may lead to general insights 
into the behavior of the explosive percolation transition, which has not been explored for 
hierarchical networks. 

II. THE BRANCHING HIERARCHICAL NETWORK AND ITS VARIANTS 

We first discuss the branching, hierarchical networks studied in this paper. The original 
network was set up as a load bearing network, and connectivity strategies were devised to 
enhance its load bearing properties [28]. Here, we describe the original 2D lattice network 
model, a critical case of the original lattice, the V— lattice [29], and its two variant load 
enhancing strategies, strategy-I and strategy-II [28]. These are shown in the four panels (a), 

4 



(b), (c), and (d) of Fig. [T] respectively. The original lattice network model has similarities 
with diverse models of real-world networks like the river network |3T] model, the voter model 
[52] . the cellular-automata model [33], and the lung inflation model [MIES]- 

The original lattice shown in panel [l|a) is, in fact, a critical case of the Coppersmith 
model for granular media [20], and also a model for the river network [21]. The lattice is a 
regular triangular lattice of sites with each site in the top layer having unit capacity. A site 
i at some layer M is connected to one of its two nearest neighbor sites chosen at random in 
the layer M + 1. The total capacity of site i in the layer M is transmitted to its left neighbor 
ii, or right neighbor v, in the layer M + 1. Thus, the capacity of the site in the layer M at 
the ith site, w{i'^^), satisfies the stochastic equation 

w{z^') = z^')w{zf-') + z'')w{tf-') + 1 (1) 

where the quantity takes the value 1 if a connection exists between «/^~^ and 

and if otherwise. The network consists of many clusters, where a cluster is defined as a set 
of sites connected with each other. The largest cluster in the system is called the maximal 
cluster, and the backbone of the maximal cluster is called the trunk. 

The original lattice has a very special realization which has the maximum trunk strength 
or capacity compared to all other possible realizations. This lattice bears a unique V— 
shaped cluster that includes all the sites in the first layer, and (M — / + 1) sites in the 
/— th layer, where M is the total number of layers. One of the arms of the V constitute 
the trunk, and all other connections run parallel to the arm of the V that is opposite to 
the trunk. Thus, this cluster includes the largest number of sites, and is thus the largest 
possible cluster the original lattice could have. We call this lattice the V— lattice (panel 
[T]^b)), and the cluster, the V— cluster. The other two variants of the original lattice which 
we have studied for simulation are the strategy-I (panel ^c)) and the strategy-II (panel 
[T]^d)) lattices. These two strategies [28] were developed for strengthening the load bearing 
capacity of the lattice by reconnecting a site on a given layer to a suitable site in the layer 
below with the reconnections restricted to at most one per layer. In strategy-I, the maximum 
number of sites are included in the maximal cluster by connecting as many disjoint clusters 
as possible to the sites on the trunk T of the maximal cluster. This is a bottom to top 
strategy. Another strategy is a top to bottom strategy where the reconnections start from 
the first layer onwards is also possible, so that the capacities of the site on the trunk in the 
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FIG. 2. (Color online) The distributions of capacities w displaying power law behavior 
P{w) = cw~°' for (a) the original lattice with exponent a = 1.3335 and = 9.3924, (b) the 
strategy-I with exponent a = 1.336 and = 3.3028, and (c) the strategy-II with exponent 
a = 1.2309 and = 0.0924. The capacity distribution for (d) the V— lattice does not display 
power law behavior. 

layer below and its subsequent sites get enhanced in every layer. This is strategy -II (see 
[25] for further details). 

We have carried out simulations for the capacity distribution P{w) of all the four lattice 
networks mentioned above. Here the capacity distribution is, in fact, the probability that 
a site has capacity w over the entire network. The distributions of capacities for the orig- 
inal, the strategy-I and the strategy-II network are governed by a power law of the form 
P{w) = cw'", where c is a proportionality constant. However, the capacity distribution 
for the V - lattice does not behave in accordance with the power law. The exponents of 
the power law for the original, the strategy-I, the strategy-II lattices are found to be 1.333, 
1.336, and 1.23 respectively as shown in the panels [2]^a), ^h), and |2](c) of Fig. |2} The 
mean-field analysis for the capacity distributions for these three networks has been carried 
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out by Kachhvah et al. [29] using the method of Takayasu [5] for the original lattice. 



III. STATISTICS OF PACKET TRANSMISSION ON THE NETWORKS 

We now consider packet transport on the hierarchical lattices. A packet or particle 
following the topology of the network hops from one site to a neighboring site along the 
link between them. In this scenario, the time the packet spends in the given site is directly 
proportional to the capacity of that site i.e. the probability Pi of a packet to be found in a 
site i of branching hierarchical network, is proportional to the capacity Wi of that site, i. e. 

Pi = ^ , (2) 

where the sum is over all the nodes of the hierarchical network. 

In statistical physics, we know that the occupation probability of a site i with energy Si 
is proportional to the Boltzmann factor 

P^ = Ce-^^% (3) 

where C is a normalization constant, /3 plays the role of the inverse temperature, /3 = 1/T, 
and Si is the energy associated with the site i with capacity Wi in the network. Equating 
the Pi in Eqs. (1) and (2), the relation between the energy associated with a site and its 
capacity Wi is given by |17| 

e-P'^ = Cw,, (4) 

If Hi be the mean occupation number of node i with its capacity Wi and where Ui < Wi, 
the total number of packets or particles |19] present on the network would then be 



N = J2{n,). (5) 



The total energy E of the network is then given hj E = Eifii. Since the packets hopping 
on the sites in the hierarchical network are distinguishable from each other, the canonical 
partition function Zj^i [IS] when distinguishable packets are present in the system is then 
given by 

Zn = J:;^^ ^e-^S^%-, (6) 

ni!n2! ■ ■ ■ Um^l 



where Yl' denotes that the indices rij in the sum satisfies the constraint = The 
factor ni\n2\ ■ ■ ■ fim^l inside the sum arises from the fact that the packets are distinguishable 
from each other. Our interest here lies only in calculating average quantities like (n^) the 
mean occupation number . These can be derived using either Zj^, or a normalized 
partition function defined as 



„2 



^ iV! = nM. ■ ■ ■ n^.V 

{■"■ii 

The grand canonical partition function [49j, using Z'j^, is given by 

oo 

E=J2e(^^^Z'^. (8) 

Af=0 

For a large number of packets, the grand-canonical partition function using either Z'j^ or 
Zn will give the same predictions for average quantities like (n^). After using Eq. (7) and 
the constraint A^ = Yli (^«) Eq. (8), the expression for the grand-canonical partition 
function, after the usual algebra, can be written as 

i=l \ni=l * / 

where the capacity Wi for the original lattice is given by the evolution equation (1). From 
the grand-canonical partition function, the mean occupation numbers [19] of packets in a 
site i of the network is calculated by using the relation 

Using equation (9) for H in the above equation, the following expression is obtained for mean 
number of packets in the site i 

S-^Wi 1 -0(ei-n)n 

(n-) = ^"=^ fll) 

After substituting the value of Ei from Eq. (3) in the above expression, the mean number 
of packets in a site i as a function of its capacity Wi is, 

where the parameter A = ce^^ is a function of the chemical potential yU. The index i here 
runs over the distinct capacities. In the limit of (nj) — )■ Wj, the nodes with high capacity 



tend to become fully occupied. For large networks, and for large capacities Wi — )■ oo, the 
above expression takes the form 



n.. 



^ ' ^ n=0 ' 



which simplifies in the large capacity limit to be 

(rii) ~ Awi. (14) 

Further from Eq. (5) for the total number of packets present in the network, and Eq. 
(13), an expression using which the chemical potential ^ oi A can be determined, is obtained 
as follows 

here the sum is over all the sites in the network. If p = N/M'^ be the occupation density, 
and Pw the capacity distribution of the network, a relation to determine the parameter A 
from Eq. (15) (after taking the sum over capacities instead of over sites) can be written as 
follows 

V P ^-^(^^^^-^" = p (16) 



w 



Eqs. (13), and (14) predict the mean occupation number {|nj) in the branching hierarchical 
networks. Since the capacities of sites in the tree-like hierarchical 2D network tend to larger 
values as the depth of network increases, taking the limit of Eq. (16) for large network size 
and large capacities, the above expression takes the form 

J2pUAw.)-p- (17) 

w 

For a network where the capacity distribution P^ is a power law i.e. P^ = cw~°', where a is 
the exponent of the power law, Eq. (17) can be used to obtain the value of A to be: 

A~(2-a)^. (18) 

We note that the original lattice, and the strategy-I and strategy-II versions discussed in 
section II have capacity distributions of the power law form [29]. 
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FIG. 3. (Color online) The mean occupation number (n^,) of sites with capacity w of (a) the 
original lattice, (b) the strategy-I lattice, (c) the strategy-II lattice, and (d) the V— lattice networks 
of (100 X 100) corresponding to the packet density p = 3. 



A. Simulation results 

We now simulate the transmission of information packets on the network. This trans- 
mission of information is modeled by packets of information hopping on the sites of the 
branching, hierarchical, heterogeneous, networks. These packets could be the data packets 
of information in the Internet, which links computers of heterogeneous and high capacities 
capable of transmitting packets at high rates. For our networks, the packets are deposited at 
a randomly chosen site on the topmost layer of the network. Each site retains the number of 
packets which saturates its capacity and the remaining packets are transmitted further. A 
packet at a given site sees the nearest neighbor sites linked to itself. If the targeted neighbor- 
ing site is not fully occupied (i.e. it has not saturated its capacity), the packet moves there 
and looks for the next site that has spare capacity. If the target site is fully occupied then 
the packet stops on the site which it occupies, and the transmission of the packet ends at 
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that site. In this fashion, all packets hop from one site to another according to the vacancy 
available on the neighboring sites, and this process continues till all the packets come to rest 
at some site. One realization of packet transmission ends here. The same number of packets 
is again placed on a randomly chosen site in the first layer of the hierarchical network, and 
the second realization of packet transmission proceeds. This stochastic process of packet 
transmission is repeated for a large number of realizations. Since our analytical expression 
for mean occupation number has the property that two nodes with the same capacity have 
the same occupation number, the mean occupation number is computed for each capacity 
instead of each site, by averaging over all the realizations of packet transmission. We finally 
look at the average occupation number (ra,) of the sites with the same capacity on the lat- 
tice, averaged over all the realizations, and compare it with the analytic result obtained in 
the earlier part of this section. 

It can be seen from Fig. [3] that the theoretical results match very well with the simulation 
results. Panels |3](a), |3](b), andjsj^c) corresponding to the original lattice, the strategy-I, and 
the strategy-II lattices respectively show a very good agreement between the simulation and 
the theoretical results obtained [50j. However, the analytic form for the capacity distribution 
for the V— lattice is unknown, so the theoretical match to it could not be shown in panel 



statistical dynamics of sites on the branching hierarchical networks. This match between 
the theoretical and simulation predictions for the original lattice network and its different 
network topologies supports our finding that the dynamics of the transmission or flow of 
the packets of information or particles in weight bearing hierarchical networks follow the 
Maxwell-Boltzmann statistics. 

IV. PERCOLATION ON THE NETWORK 

We study site percolation on the hierarchical networks set up here due to its relevance 
for the problem of information transfer on the lattice. The site percolation problem is set 
up as follows: If at any time, information packets occupy/do not occupy, a given site in 
the network, that site is said to be occupied/unoccupied at that time. The network then 
can be thought of being comprised of two sub-lattice networks, one being the unoccupied 
or free sub-lattice network, and the other being the occupied sub-lattice network. Here, we 




Eq. (13), along with the Eq. (14) predicts the mean occupation number 
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FIG. 4. (Color online) The percolation strength as a function of packet density /i, averaged 
over 500 realizations, for the unoccupied sub-lattice in (a) the original lattice, (b) the strategy - I 
lattice, (c) the strategy - II lattice, and (d) the V— lattice networks of (100 x 100) sites. 

are interested in studying the transition to percolation in different lattice networks. From 
site percolation theory it is known that if the sub-lattice network of free or unoccupied 
sites percolates, it implies that a finite fraction of remaining unoccupied sites form a single 
connected giant cluster. We can identify this the state in which the single giant 

cluster could make communication possible throughout the network. If the free sub-lattice 
network does not percolate, the free or unoccupied sites then form mutually disconnected, 
and fragmented small clusters instead of a single giant cluster. In this state, communication 
through the network is not possible. It is interesting to study the transition to percolation 
in the unoccupied or occupied sub-lattice network. We anticipate a transition from the 
percolating to the non-percolating state when the packet density /i changes. In order to 
analyze the transition from the percolating to the non-percolating state in networks, we 
simulate the stochastic dynamics of packet transmission in the original lattice, the strategy- 
I and the strategy-II lattices, and the V— lattice for different values of the packet density 
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FIG. 5. (Color online) The susceptibility x as a function of packet density fi, averaged over 500 
realizations, for the unoccupied sub-lattice in (a) the original lattice, (b) the strategy - I lattice, 
(c) and the strategy - II lattice networks of (100 x 100) sites. 



yU. Two types of transition are seen in the system, as we will see in the next section. 



A. The two transitions: Percolation and explosive percolation 

To study the phase transition, we define the order parameter to be the percolation 
strength S = Sm/L [38], where Sm indicates the number of sites belonging to the largest 
connected cluster, and L = M^, the total number of sites in the network. We computed the 
percolation strength 5* of the largest cluster of the free sub-lattice network as a function of 
the packet density fi averaged over a number of realizations as shown in Fig. |4j The panels 
|4|^a), |4](b), |4]^c), and |4]^d) display the percolation strength as a function of the packet 
density /i for the original, the strategy-I, the strategy-II lattices, and the V-lattice respec- 
tively. The first three panels |4]^a), ^h), ^c) corresponding to the original, the strategy-I, 
and the strategy-II lattice network demonstrate that the transition from percolating to a 
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FIG. 6. (Color online) The percolation strength 5 as a function of packet density displays 
power law behavior S oc (/i — /^c)"'^ corresponding to 500 realizations, for the unoccupied sub- 
lattice in (a) the original lattice with rj = 0.0466, = 0.00004, (b) the strategy - I with i] = 0.158, 
= 0.00108, and (c) the strategy - II lattice with r] = 0.368, = 0.00143, networks of L = 10^ 
sites. 

non-percolating state in these lattices is actually a continuous phase transition. In each 
panel of Fig. |4| fic is the critical packet density where a transition from a percolating to a 
non-percolating state occurs. 

The most striking result is however observed for the critical case of the original lattice 
the V— lattice [29], in which a sharp, discontinuous percolation transition appears unlike 
the other three cases. 

B. Finite size scaling and critical exponents: 

The percolation strength S of the largest cluster of unoccupied sites when plotted versus 
(/i — fic), shows that 5* displays power law behavior of the form S* ~ (/i — fic)^^ for the 
all the three lattice network indicating the existence of a continuous percolation transition 
as shown in Fig. [6} except for the case of the critical V— lattice. For the V— lattice, the 
percolation transition is discontinuous, and hence scaling does not work. The exponents rj 
of the power law for the original, the strategy-I, and the strategy-II lattices are found to be 
0.466, 0.158, and 0.368 respectively. 

For continuous percolation transitions, the percolation strength S* of a network of L sites, 
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FIG. 7. (Color online) The percolation strength S at the critical point /ic, as a function of number 
of sites L and corresponding to 500 realizations, displays power law behavior 5 ~ L^/^ ^ for the 
unoccupied sub-lattice in (a) the original lattice with fi/v = 0.7902 it 0.0023, (b) the strategy - I 
with P/v = 0.9859 ± 0.015, and (c) the strategy - II networks with 13 /v = 0.9913 ± 0.0041. 
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FIG. 8. (Color online) The susceptibility x the critical point function of number of 

sites L and corresponding to 500 realizations, displays power law behavior S ~ L~^l^ ^ for the 
unoccupied sub-lattice in (a) the original lattice with ^jv = 0.31 ±0.0271, (b) the strategy - I with 
-f/u = 0.1063 ± 0.0092, and (c) the strategy - II networks with -//u = 0.055 ± 0.0161. 

obeys the following finite size scaling relation 

S = L^/''f[{f,-fx,)L'/% (19) 

where P is the critical exponent for the order parameter, u that of the correlation length, and 
/(■) is a universal function. We also define the susceptibility x = \\j (S^) ~ which 
quantifies the amplitude of the fluctuations of the percolation strength S. The susceptibility 
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FIG. 9. (Color online) The log - log plot of critical packet density fic as a function of the network 
of L sites for (a) the original lattice with — = 0.36 it 0.13, (b) the strategy-I lattice with 
= 0.43 ± 0.02, (c) the strategy-II lattice with -l/u = 0.437 ± 0.026. 



in our cases, shows finite size scaling, with the scahng function 

X = L~'/''9[{^^-^^e)L'/% (20) 

where 7 is a critical exponent of the susceptibility and g{-) a universal function. The 
susceptibility x is used to determine the critical point fic- The panels |5](a), [Sj^b), and|5](c) 
display the susceptibility x against the packet density fi for the original, the strategy-I, 
and the strategy-II lattice respectively. The critical point fic in these cases is determined 
from the first minimum after the peak (for the unoccupied sub-lattice networks) as shown 
in panels [5|a), |5](b), and|5](c). 

It is known that [38j, the susceptibility, and the percolation strength at critical point 
are directly related, and their critical exponents satisfy the equality P/u + 'y/u = 1. For 
the original, the strategy-I and the strategy-II lattices, where a continuous phase transition 
is observed, we computed the critical exponents by examining the finite size scaling of the 
two variables, the susceptibility x and the percolation strength S at the critical point and 



comparing these with Eqs. (19) and (20) (see the panels of Fig. 10 and Fig. |8]). From the 
panels [Io|^a) and [sj^a) for the original lattice, it can be seen that the values of exponents 
P/u and 7/1/ fall between and 1, and their sum P/u + 7/z/ is very close to 1. The same 
can be observed for the strategy-I from the panels flOtb) andlstb), and for the strategy-II 



lattice from the panels 10 c) andpTc). If the transition is discontinuous as in the case of 



the V- lattice, the finite size scaling law from Eq. (19) does not hold, and the percolation 
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transition is of the explosive type. Since for the V-lattice, 5* at critical point does not vanish 
in the infinite size limit, the ratio = 0, and the ratio ^jv = 1 as the susceptibility x is 
extensive (from panel [sj^d)) [38] . Thus the relation between the critical exponents is satisfied 
for this case as well. 

It is also observed that the critical packet density fic{L) corresponding to a network of L 
sites obeys the relation 

/i,(L) =/ie + L-l/^ (21) 



as can be verified from the panels|9|^a),|9|^b),[9|^c) for the original, the strategy-I, and strategy- 
II lattice networks, respectively. This expression can be used to compute the values of u 
corresponding to the three lattices. The values of all three exponents (3, 7 and u are tabulated 
in Table I for the original, the strategy-I and the strategy-II lattices. 

Table I: The values of the critical exponents /3, 7 and u of the percolation transition for 
the original, the strategy-I and the strategy-II lattices. 



Network 


/3 


7 


1/u 


/3/z/±7/z/^ 1 


Original 


2.19 ±0.0064 


0.86 ±0.0753 


0.36 ±0.13 


1.1002 


Strategy-I 


2.27 ±0.0348 


0.256 ±0.0214 


0.43 ±0.02 


1.0922 


Strategy-II 


2.265 ± 0.0094 


0.114 ±0.0368 


0.437 ±0.026 


1.0463 



The scaling behavior of the percolation strength S for the original, the strategy-I, and the 



strategy-II lattice networks are shown in Fig. 10 



We, further computed the size of the largest cluster Sm at the onset of the phase transition 
i. e., before the giant cluster is formed (at the jump in case of the V-lattice), as a function 



of the network size. From the three panels of Fig. 11 corresponding to three different 
topologies, the size of the largest cluster Sm prior to the transition threshold does not grow 
proportionally to the network size M, instead it obeys a power law of the form of Sm ~ M'^. 
We note that in the case of the original lattice ( = 1.6, whereas for strategies I and 11, ( ^ 2 
indicating that the strategies have been successful in incorporating a substantial fraction of 
the sites of the lattice in the maximal cluster. 
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FIG. 10. (Color online) Scaling of percolation strength S corresponding to different network sizes 
M i. e., S/M^ versus fj,/M^ for (a) the original lattice with x = 2.7, y = 1.62, (b) the strategy-I 
lattice with x = 2.7, y = 2, (c) the strategy-II lattice with x = 3.1, y = 2. 




FIG. 11. (Color online) The log - log plot of Sm, the largest cluster size at the onset of percolation 
transition (at the jump in case of the V-lattice), as a function of the network size M for (a) the 
original, (b) the strategy-I, (c) the strategy-II, and (d) the V— lattice networks. 
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V. CONCLUSIONS 



The study of the transport of packets on hierarchical networks has yielded some intriguing 
insights into the statistics and the percolation behavior of such systems. Most earlier studies 
of transport on networks have been performed on either scale-free or random networks and 
not on hierarchical ones. The models used by us to study traffic of packets, are typical 
example of tree- or river- like uni-directional branching hierarchical networks, as in computer 
networks with backbones in which the capacity distribution of sites is more significant than 
the degree distribution for the analysis of traffic. These networks are highly heterogeneous, 
and display power law behavior for the capacity distribution. If distinguishable packets are 
transported on these networks, an analytical expression for the mean occupation number of 
sites on a hierarchical network can be obtained as a function of its capacity distribution, 
which matches with the numerical simulations, and confirms that the packet distribution 
satisfies the Maxwell-Boltzmann distribution. 

Site percolation studies on networks of a given geometry and capacity distribution can also 
yield important insights into their efficiency for information transport. Hence, we studied 
the site percolation problem for different topologies of the hierarchical network. We have 
seen that the percolation transition in the original model and its two modified versions is 
a continuous transition, whereas the percolation transition is explosive (discontinuous) in 
the critical case of the original lattice, i.e. the V— lattice. We have recently seen that 
the avalanche distribution on the V— lattice shows power law behavior, unlike the peaked 
behavior seen on the other lattices. Thus the V— lattice is a special lattice, which promotes 
rapid transport of traffic. Recently, the explosive transition has been explored extensively in 
various scale-free and random networks. Some results indicate that the explosive percolation 
transition, for certain classes of networks, is not a genuine discontinuous transition, and the 
discontinuity is a consequence of the finite size of the lattice. However, these lattices are quite 
different from our hierarchical lattices, and the question of whether the explosive percolation 
on the V— lattice will survive the infinite size limit is worth independent study. We hope 
to explore this question, as well as the question of the optimality of the V— geometry for 
efficient transport in future work. 
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